Crystal nucleation and cluster-growth kinetics 
in a model glass under shear 
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Crystal nucleation and growth processes induced by an externally applied shear strain in a model 
metallic glass are studied by means of nonequilibrium molecular dynamics simulations, in a range 
of temperatures. We observe that the nucleation-growth process takes place after a transient, in- 
duction regime. The critical cluster size and the lag-time associated with this induction period 
are determined from a mean-first passage time analysis. The laws that describe the cluster growth 
process are studied as a function of temperature and strain rate. A theoretical model for crystal- 
lization kinetics that includes the time-dependence for nucleation and cluster growth is developed 
within the framework of the Kolmogorov-Johnson-Mehl-Avrami scenario and is compared with the 
molecular dynamics data. Scalings for the cluster growth laws and for the crystallization kinetics 
are also proposed and tested. The observed nucleation rates are found to display a nonmonotonic 
strain rate dependency. 

PACS numbers: 46.35.+Z, 05.70.Ln, 64.60.qe, 64.70.pe 



I. INTRODUCTION 

The study of phase transformation between liquid 
and crystal through a nucleation and subsequent growth 
regime is a problem with a very long history [l], Q ■ Un- 
derstanding the crystallization process, including the rate 
of phase transition and the morphology of the crystal 
formed, has a great importance for many technological 
applications. The situation becomes even more complex, 
when the phase transformation takes place in a system 
under external drive due to shear flow, electric, magnetic 
or laser-optical fields etc. (for a recent review see (H). 
The problem belongs then to the class of nonequilibrium 
processes in driven materials, which has attracted atten- 
tion more recently 

Numerical simulations techniques have allowed one to 
obtain a series of important results for the kinetics of 
crystallization in systems driven by an imposed shear 
now , including glasses [10- 12], semicrystalline poly- 
mers [13|, colloidal suspensions [l4 ITo| and Ising model 
[l^ . The generally established outcome here, which con- 
firms the experimental observations (see, e.g., [l7l - l20l |). is 
that the shear drive can have a significant impact on the 
various aspects of the nonequilibrium phase transitions, 
in particular, on the transition, nucleation and crystal 
growth rates as well as on the induction time 0,[H|. 

In terms of nucleation, the influence of a finite shear 
rate on the structural ordering of a system appears to 
be, in general, that a small shear rate speeds u p n ucle- 
ation, while larger shear rates prevent ordering [lot [l6| . 
Fluid-crystal coexistence can also be affected by shear, 
as found in Ref. [3] where crystallization is shown to be 
suppressed by flow. Hence, the shear flow influences both 
the thermodynamic and the kinetic aspects of nucleation, 
in a way that may depend on the depth of supercool- 
ing and on the intensity of the strain rate, presumably 



compared to the system internal relaxation time. For a 
deeply supercooled, glassy system, the relaxation time is 
essentially infinite, so that a finite shear rate will always 
have a favorable impact on nucleation. 



The cluster growth process, which follows nucleation, 
is also expected to be affected by a finite strain rate. For a 
specific, Ising like, two-dimensional system, it was shown 
recently that a moderate shear-drive plays a significant 
role in crystal erosion and growth governing b y s ingle 
"particle" attachment and coalescence processes [15|, [l6| ■ 
More generally, one may inquire how the domain growth 
law are affected by the shear rate. In systems at rest, the 
standard description for the apparition of a crystalline 
phase through nucleation and cluster growth, is generally 
reproduced by the classical Kolmogorov-Johnson-Mehl- 
Avrami (KJMA) theory [l]. In the present work, we ap- 
ply the extension of the KJMA theory for the case of the 
time-dependent nucleation and growth in a model metal- 
lic glass under shear drive. Extension of the theory is 
tested together with nonequilibrium molecular dynam- 
ics simulation data at different temperatures and wide 
range of shear rates. The data is compatible with the 
extension of the theory that includes a finite lag-time for 
appearance before the onset of steady state nucleation. 



The numerical model used in our simulations, as well 
as the extended KJMA theory, are presented in sections 
ITTland llTTl The results of molecular dynamics simulations 
and parameters of the KJMA theory, in particular, the 
critical cluster size, the lag-time and the steady state 
nucleation rate, are analyzed and discussed in section 
HVl Finally, we present our conclusions in section [V] 
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II. MODEL, SIMULATION DETAILS AND 
CLUSTER ANALYSIS 



In this work we focus on a system of particles interact- 
ing through the spherically symmetric Dzugutov poten- 
tial 



U(r*)/e = A (r*- m -B) exp 



H{a-r*) 



B exp 
r/a, 



-b 



H(b-r* 



(1) 



where a and e define the unit length and energy, re- 
spectively 1 . The parameters A, B, m, a, b and c are 
chosen as proposed originally in Ref. [22|], the Heav- 
iside step function H(. . .) sets the range of the vari- 
ous contributions in relation (JTJ) . Besides a minimum 
LT( r min = 1.13) = — 0.581e, the potential includes a max- 
imum £/(r* nax = 1.628) = 0.46e and falls off rapidly with 
the interatomic distance r. Such a short-ranged, oscillat- 
ing interaction mimics in a simple way the ion-ion inter- 
action influenced by the electron screening effects in the 
metallic systems. Moreover, the maximum in potential 
([1]) reflects the first of the Friedel oscillations and favors 
icosahedral local order in the system, therefore making 
it a good glass- former at low pressure [H| . 

The system under study and the simulation setup are 
completely identical to the considered one in Ref. fl2^ . 
Namely, the system consists of TV = 19652 particles 
within the simulation box L 3 with L = 28.55cr that 
corresponds to the density p = 0.84er~ 3 . A set of 
glassy samples is prepared by fast cooling from the equi- 
librated liquid state to the temperatures T — 0.01, 
0.03 and 0.06e/fcB that is well below the melting point 
T = 0.5e/k B 0. 

The shear drive is applied by moving two amorphous 
walls created at the sides of the simulation cell perpen- 
dicular to the e y direction. The bottom wall is fixed, 
whereas the top wall is moving in the x-direction with 
the instantaneous velocity u(i) = yLtyex at a constant 
strain rate 7 and pressure P yy — 7.62e/er 3 , which in the 
equilibrium phase diagram would favor the fee phase. 
Here L(t) is the instantaneous distance between the walls. 

To identify the nuclei of the ordered phase (clusters) 
we use a cluster analysis, which is based on the consid- 
eration of the local environment around each particle by 
means of a (2 x I + l)-dimensional complex vector with 



1 For convenience, all quantities are expressed in reduced form. 
The time unit is r = Uv/mo/e, mo is a particle mass, the strain 
rate is in units of r , the temperature T is in units of e/ks 
and the pressure is in units of e/cr 3 , where ks is the Boltzmann 
constant. 



the components 



qim(i) 



N b (i) 



(2) 



where Yi m (9ij,<fij) are the spherical harmonics with the 
polar 8ij and azimuthal tpy angles between radius- vector 
Yij and a reference direction; iV&(£) is the number of 
neighbors for a particle i, which are the particles located 
within a sphere of the radius |r.y| = 1.5a (see Ref. [HI). 
Following the ten Wolde-Frenkel scheme [26|], we specify 
the pair of neighbors, particles i and j, as correlated into 
an ordered structure if the following condition is satisfied: 



96m(i)<76m(j) 



> 0.5, 



where the normalization 

qim(i) = 



qim{i) 



1/2 



(3) 



(4) 



sets the maximum possible value in the r.h.s. of inequal- 
ity © equal to unity. Moreover, to exclude from consid- 
eration the structures with a negligible number of bonds 
per atom, which occurs even in liquid phase, we apply 
the following additional restriction [26[ : particle i is con- 
sidered as included into a crystalline structure if it has 
seven and more neighbors satisfying the condition (j3|). 



III. NUCLEATION AND GROWTH KINETICS 

According to the KJMA theory for crystallization ki- 
netics |l| , the fraction of material transformed into a crys- 
talline phase at a given time t is defined by 



a(t) = 1 — exp 







,(t',t)dt' 



(5) 



where I(t) is the nucleation rate and v ex (t',t) is the vol- 
ume at time t of a nucleus formed at time t': 



V ex (t',t) = Cg 



G(t")dt" 



(G) 



G{t) is the growth rate of the nucleus radius, c g is a di- 
mensionless shape factor. This description is obviously 
correct if critical sized nuclei grow isotropically and are 
much smaller than the system size. In the simplest ver- 
sion of the theory, the growth rate does not depend on 
size or time, G(t) ~ Q c and the nucleation rate I(t) is 
approximated by the steady-state nucleation rate I s , - 
Eqs. (O and ^ are simplified to give a well-known ex- 
pression for steady-state homogeneous nucleation kinet- 



a(t) = 1 — exp — 



(7) 
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However, if the time scale of the transient regime, which 
precedes the steady nucleation and growth kinetics, is 
comparable to nucleation and growth time scales, a lag- 
time t c that accounts for the non-stationary character of 
the transition kinetics should be introduced, as discussed 
below. 

The growth law of a supercritical cluster is commonly 
chosen (see Ref. [if, p. 378) to be of the form 

R(t) = (g c ty, (8) 

with the growth rate averaged over directions 2 

G(t) = vG" c t v - x . (9) 

Here, R is the averaged radius of a crystallite, the growth 
constant Q c and the exponent v take positive values and 
are determined by the growth mechanisms. Note that 
the term Q c has a dimension of (length) ^ (time) -1 . 

Taking into account the last equation, one obtains the 
volume of the supercritical cluster V(t) — c g R(t) 3 = 
Cg(G c t) 3 "- As a result, the growth law of a supercriti- 
cal cluster can be written as 



N{t c ,t) = N c + c gPs g?{t-t c ) 



\3i/ 



t > tr 



(10) 



where N c is a critical cluster size, p s is a numerical den- 
sity of the crystalline cluster and t c is the mean lag-time 
for the appearance of a critical cluster. This equation 
can be rewritten in the dimcnsionless form: 



, _ c g p s (g c t c ) 3L ' ^ . 3v 



(11) 



where 



£ = -, £>1. 

Assuming that Eq. (fTU|) holds to describe the growth of 
a supercritical cluster, we obtain equation for the time- 
dependent extended volume of a single cluster 



Vex (j'ct ^) 



N{t c ,t) 



Ps 



(12) 



1 



-[N c + CgPs g 3 c v (t-t c ) 3v ] 

Ps 



Further, the simplest model to take into account the 
existence of a transient regime on the nucleation rate 
consists in assuming that this rate is zero until t Cl and 
becomes constant hereafter [H, [29| . This corresponds to 
a function I(t) given by 



I(t) = I s H(t-t c ), 



(13) 



2 Generally, clusters can have a distribution of shapes and struc- 
tures. Here we use a simplified description in terms of appropri- 
ately averaged cluster shapes [27l . 



where H(t) is the Heaviside step function. By inserting 
Eqs. (fT2"j) and ([TBI into Eq. ([5]), we obtain for nucleation- 
growth regime t > t c the following equation: 



a(t) 



= 1 — exp 



I .Nr. 



+ 



(t - t c ) 

{3v + l)N c 



(14) 



(t-t c ) 



3v+l 



where the normalization factor indicates the possi- 
bility of incomplete crystallization of the parent phase, 
< < 1. By analogy with the growth law for a su- 
percritical cluster [see Eq. (jlll) ] the last equation can be 
also written in the dimcnsionless form: 



«(0 



1 — exp 



(e-i) 

{Zv + l)N e 



(15) 



(£-1) 



3H-1 



When the critical cluster size N c is much smaller than 
the system size, the first (linear) term in the exponent of 
Eqs. (|14p and (fl~5|) can be neglected. Moreover, if the lag- 
time t c is negligible in the nucleation and growth kinetics, 
then Eq. (|14|) reduces to the well-known Avrami equation 
[Eq. ©]. 



IV. RESULTS 

A. Critical cluster 

The critical cluster size iV c is one of the crucial param- 
eters in nucleation theories. Since the critical clusters 
are undetectable by the common traditional experimen- 
tal tools, especially at the high supercoolings correspond- 
ing to a glassy phase, advanced methods must be used to 
define the critical cluster size and to clarify the question 
about subcritical cluster morphology (see review [30|). 
Here one can mention the advanced Koster's method [31j 
and the way based on the accurate study of static struc- 
ture factor data 



321 ] . On the other hand, molecular dy- 



namics simulations allow one to identify clusters of all 
sizes including supercritical solid clusters of the nucle- 
ation regime as well as subcritical clusters of the tran- 
sient regime. Therefore, the critical cluster size can be 
identified, if a method to define correctly the boundary 
between these regimes is found. 

Yasuoko-Matsumoto method. - A first method for 
identifying N c is based on the consideration of the time- 
dependent total number of clusters whose size is larger 
than a given value N* , i.e. 



AW(t) 
f(N*,t)= £ n s (t), 



(16) 



=N* 



where n s (t) is the time-dependent cluster size distribu- 
tion, A r max (t) = max Si „ s ^o[ n s(0] i s l ne srze 01 ^ ne largest 
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Figure 1: (Color online) Main: Time evolution of the total 
number of clusters f(N*,t) with sizes larger than the thresh- 
old value N* . The system is under shear with the strain rate 
7 = O.OOOlr -1 at the temperature T = O.Ole/fcs. Here N* 
increases by steps of 3 from the top curve with iV* = 3 to 
the bottom curve with N* — 47. Inset: Full circles present 
the slope of the curves f(N*,t) at the different values of N* 
in the time range indicated by the vertical lines in the main 
figure, corresponding to nucleation regime. The solid line is 
the fit by the function s(N*) = so + aexp(— bN*) with the 
parameters sq = 0.01, a = 0.15 and 6 = 0.06 (see Fig. 2 ir 
Ref. HH] for comparison). 



cluster at the time t. Obviously, at N* = N c Eq. flgf 
defines the total number of supercritical clusters f(N c , t) 
formed at time t, whereas (l/V)df(N c , t)/dt is the nucle- 
ation rate V is the volume. In the steady nucleatior 
regime, the rate I(t) = const = I s , defined by the slope oi 
f(N c ,t)/V, is independent of N* for the range of super- 
critical clusters, i.e. N* > N c . Note that this is correct 
only if the cluster growth rate is independent of time anc 
size. As a result, the time-dependent curves f(N*,t) al 
different N* > N c must be simply shifted and have the 
same slope for the steady nucleation regime. This reg- 
ularity will appear in the vicinity to the critical cluster 
size N c [H. 

To test the suitability of this method for extraction of 
the critical cluster size N c , we compute f(N*,t). The 
time evolution of f(N*,t) at the different threshold val- 
ues of ./V* is shown in Fig. [T] As expected, the curves 
shift in t with the increase of N*. A linear growth in 
the nucleation regime (marked one in the main figure) 
is observed for all the cases including the case with ex- 
tremely small threshold value N* — 3. At the same time, 
the expected regularity that would indicate the indepen- 
dence of the f(N*, i)-slope on N* is never observed (see 
inset of Fig. [1]). Instead, we find that the dependence of 
the slope on N* is well fitted by an exponential decay 
[34| . Therefore, we conclude that in the present case, the 
method does not allow one to find explicitly the values 
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Figure 2: (Color online) Mean first-passage time distributions 
multiplied by the corresponding strain rates. The system is at 
the temperature T = 0.01e/fcs. (a) Circles correspond to MD 
simulations data as averaged over set of 50 independent runs 
at the strain rate j = 0.008 and O.Olr - . Solid lines are their 
fit to Eq. (|18[) . Dashed lines mark for both cases the criticai 
cfuster size defined through the iocation of the inflection point 
of the fitting curves, (b) MD simuiations data at the different 
strain rates j G [0.000 1, 0.01]r _1 . 
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Figure 3: (Color online) Growth of the largest cluster with the 
strain 7 = -yt at the different strain rate 7 G [0.000 1, 0.01]r _1 
and the temperature T = O.Ole/fcs. The strain rate grows 
from left to right for a plot. Dashed line corresponds to the 
value 20. Full circles indicate the critical cluster sizes defined 
as described in the text. 



of the critical cluster size N c . On the other hand, such 
a behavior of f(N*,t) can arise because both nucleation 
and growth take the comparable time scales [HI, and 
the growth rate is a time- and size-dependent. In that 
case, a supercritical cluster is growing by a single-particle 
attachment as well as by a cluster coalescence, where a 
larger cluster merges with a smaller one. 
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Figure 4: (Color online) Cluster growth curves of a sheared system at the temperatures T = 0.01, 0.03 and 0.06e/fcs and the 
different strain rates (solid lines) rescaled according to Eq. (|lip onto the critical cluster size N c and the lag-time t c . Full dots 
present the parabolic fit A(£, — l) 2 , where A = 1.7 ± 0.3. 



Mean first-passage time method. - This method fo- 
cuses on the evolution of the largest cluster with the aim 
to define the average time of first appearance of a cluster 
with size N (see Refs. (HHl]), i.e. 

i M 

»=i 

where rjy is the time of the first appearance of the N- 
sized cluster after a single run and M is the total number 
of runs, i = l,2,3, ...M. If the nucleation is followed by 
fast cluster growth, the f(N) has a pronounced sigmoidal 
form and can be fitted by 

t(N) = t c [l + erf((iV - N c )c)], (18) 

where erf (. . .) is the ordinary error function, c defines 
the curvature and is related to the Zeldovich factor 
Z = cj \pK [36]. Then, the critical cluster size N c can 
be simply defined by the position of the inflection point 
in MFPT [38|, which indicates the onset of stable clus- 
ter growth, whereas the term r(JV c ) will characterize the 
mean lag-time t c for the appearance of a critical cluster. 
It is important to note that if the transient regime is in- 
significant and the transition is characterized by steady- 
state nucleation mainly, then MFPT method allows one 
to estimate the nucleation rates directly as an inverse 
height of plateau in MFPT divided by volume [13, S^. 

Figure [5] shows MFPT distributions averaged over set 
of independent runs, where each curve corresponds to 



a fixed temperature and strain rate. Although a fit of 
Eq. (ITBl to the data is suitable and the plateau in MFPT 
is observable for all cases [Fig. [2fb)], the position of the 
plateau is difficult to locate accurately [Fig.[2fa)]. Such a 
behavior indicates once again the possible non-stationary 
character of the cluster nucleation-growth process and 
that nucleation and growth occur at the comparable time 
scales j35|. Nevertheless, as can be seen from Fig. (Ha), 
the inflection point in MFPT distributions, which is as- 
sociated with the critical cluster size N c and the lag-time 
t c , is well-defined. 



B. Cluster growth 

Figure [3] shows the growth curves of the largest 
cluster in the system under shear, for strain rates 
7 e [0.000 1, 0.01]t _1 , and at a low temperature T = 
O.Ole/fcs. It is seen from the figure that all the curves 
indicate a steady growth for a cluster with size larger 
than N = 20 particles. Therefore, the threshold value 
of a cluster size, associated in CNT with a critical clus- 
ter size N c , at which the steady growth starts, must be 
relatively small. The values of iV c defined by means of 
MFPT method are also presented in Fig. [31 These small 
values of N c are qualitatively in agreement with CNT, 
which predicts a decrease of the critical cluster size with 
the supercooling. Another interesting observation is re- 
lated with the existence of a transient regime, which pre- 
cedes the nucleation and growth processes and causes 
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the delay of the cluster growth. The transien 
can be characterized by the lag-time t c , whicl 
the time required for the appearance of a criti 
cluster [H, [33[. Hence, the shift of the growtl 
observed in Fig. [3J indicates directly a 7— deper 
the lag-times. We note that the similar growth c 
also observed for the cases with T = 0.03 and 
Growth laws of the largest cluster at varioi 
rates and temperatures are presented in Fig. 21 1 
scaling form described by Eq. (dJ). First, as can 
all curves for the fixed temperature are collapse 
single master curve, which indicates the univer 
acter of the growth kinetics. Moreover, the mas 
is very well fitted by Eq. (fTT|) with the growth < 
v = 2/3 and the factor c g p s (Q c t c ) 2 /N c = 1.7 : 
all considered temperatures. Both parameters a 
be T-independent for the temperatures we have 
On the other hand, it appears that the growth 
Q c correlates with the lag-time t c and the critic; 
size N r , and so we have 
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Taking into account that v = 1 for 3D uniform crystalline 
growth controlled by interface transfer, the smaller value 
of the growth exponent, v — 2/3, can reflect the influ- 
ence of shear on the growth mechanism, where the corre- 
sponding cluster growth is considered as an averaged one 
over directions. Note that a crystal growth law oc t Zv 
with a small exponent 3i/ = 1, was observed in the diffu- 
sion wave spectroscopy "echo" experiments for colloidal 
glasses under shear of Ref. 18]. While we do not have 
any theoretical explanation for the empirical correlation 
(fTTl)) . this correlation is clearly associated with the influ- 
ence of shear on the kinetic aspect of the nucleation and 
growth process. The growth constant Q c is defined by 
a particle attachment frequency, whereas the lag-time t c 
characterizes shear-induced "unjamming" of the glassy 
system (T^ |. 



C. Crystallization kinetics and nucleation rates 

We now come to a discussion of the crystallization 
kinetics of the glassy system under strain. The time- 
dependent crystalline phase fraction ot(t), as resulting 
from cluster analysis applied to our simulation data at 
three fixed temperatures T = 0.01, 0.03 and 0.06e//cs 
and four fixed values of the strain rate 7, is presented 
in Fig. [5] The data for each case presented in the figure 
is the result of averaging over a set of independent runs. 
One can see from the figure, that the evolution of the 
crystal fraction is characterized by three distinct regimes, 
in analogy with the transition without external drive. 
In the first regime the crystalline fraction is practically 
negligible. The characteristic time scale for this stage is 
defined by the time elapsed between the system quench 



Figure 5: (Color online) Time evolution of the crystalline 
fraction at temperatures T = 0.01, 0.03 and 0.06e/fcs and for 
different strain rates. In all plots, each curve corresponds to a 
fixed value of 7 = [0.05, 0.001, 0.005, 0.0001]t-\ The strain 
rate grows from right to left for a plot. The circles are results 
obtained from simulations, solid curves present the fits with 
Eq. (|14l) . For clarity, the curves are presented only for four 
values of strain rate from the seven considered. 




Figure 6: (Color online) Time evolution of the crystalline 
fraction at a temperature T = O.Ole/fcs and strain rates 7 £ 
[0.000 1; 0.01]r _1 . Curves are rescaled according to Eq. (fT5l) . 
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Table I: Parameters of the crystallization kinetics: the system temperature T, the strain rate 7, the critical cluster size N c , the 
lag-time i c , the product of cluster growth characteristics c g p a Q^. and the crystallization factor a x . The numerical density of 
the amorphous phase is p am = 0.85<r -3 and of the crystalline fee phase is p c = 1.09<r -3 . 



T (e/k E 



7 (r- 1 ) 



tc (t) 



0.01 
0.01 
0.01 
0.01 
0.01 
0.01 
0.01 



0.000 1 
0.000 5 
0.001 
0.002 
0.005 
0.008 
0.01 



10 ±1 
11 
15 

20 ±1 

21 ±2 
19 ±1 
16 ±1 



1 400 ± 350 
310 ±10 
235 ±5 
195 ±5 
200 ±6 
215 ± 2.5 
280 ±5 



0.8 
19.4 
46.1 
84.4 
89.3 
69.7 
34.6 



0.77 ±0.01 
0.68 ±0.01 
0.68 ±0.01 
0.69 ±0.01 
0.69 ±0.01 
0.7 ±0.01 
0.7 ±0.01 



0.03 
0.03 
0.03 
0.03 
0.03 
0.03 
0.03 



0.000 1 
0.000 5 
0.001 
0.002 
0.005 
0.008 
0.01 



9 

13 ±1 
15 ± 1 
21 ±2 
20 ±2 
18 ±1 
15 ± 1 



850 ± 150 
298 ±6 
235 ±5 
195 ±5 

176 ± 10 
178 

179 ± 15 



2.1 
24.8 
46.1 
93.9 
109.6 
97.0 
79.5 



0.77 ± 0.02 
0.71 ±0.01 
0.77 ±0.01 
0.76 ±0.01 
0.75 ±0.01 
0.78 ±0.02 
0.75 ± 0.01 



0.06 
0.06 
0.06 
0.06 
0.06 
0.06 
0.06 



0.000 1 
0.000 5 
0.001 
0.002 
0.005 
0.008 
0.01 



13 ± 1 

14 ± 1 
16 ± 1 
18 ±1 

15 ± 1 
14 
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Figure 7: (Color online) Main: Nucleation rate isotherms as 
function of the strain rate 7 at the temperatures T = 0.01, 
0.03 and 0.06e/fc_g. Full symbols depict the data extracted 
from crystallization kinetics; open symbols present the results 
obtained according to Yatsuoko-Matsumoto method. The 
short solid line indicates the unit slope. Inset: Inverse lag- 
time vs. strain rate at temperatures T = 0.01, 0.03 and 
0.06e/fcs- The correspondence between curves and tempera- 
tures is the same as in the main figure. 



and the formation of a critical cluster. The second regime 
corresponds to the growth of the crystalline fraction. Fi- 
nally, in the third regime, the crystalline growth is essen- 
tially terminated, and the small increase app ears due to 
coarsening and defect removal processes Moreover, 
Fig. [5] shows the fit of these results by the theoretical 
model (ITU) presented in Sec. IIIII Only two adjustable 
parameters, and I s , are needed to reproduce simu- 
lation results for all considered cases, as the values of 
all other terms included in Eq. (|14p are known from the 
cluster growth curves discussed above. It is clear that 
the parameter aioo defines the final part of the crystal- 
lized fraction and, thereby, can be found from the final 
plateau of a(t). Then, if theoretical model (TT4|) is capa- 
ble to reproduce the data, then we obtain within such 
a model an additional tool to extract the steady-state 
nucleation rate I s from crystallization kinetics data 3 . 

As can be seen from Fig. [5J an excellent agreement be- 
tween the molecular dynamics simulation data and the 
model (TT4|) is obtained for all considered temperatures 
and values of the strain rate 7. The growth exponent 
v was taken as 2/3, the product c g p s Gci th e lag-time t c 
and the critical cluster size N c were used immediately as 
found from cluster growth analysis [see Table [TJ . More- 



3 A similar method was used in Ref. [2S| to estimate the growth 
and nucleation rates. 
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over, from Fig. [5j the values of do not saturate to 
unity. This indicates incomplete crystallization of the 
glassy system, ~ 65 4- 70% of the bulk, and represents a 
consequence of the shear drive. Note, that a similar ef- 
fect was observed by Rottler-Srolovitz for shear-induced 
alignment in polycrystalline bilayer systems, although 
the final ordered fraction was found to be less than in 
our case [39j . 

The slope of a(t) in the nucleation-growth regime ap- 
pears to be the same for all the considered shear rates 
and temperatures, indicating the universal character of 
the crystallization kinetics. In Fig. [51 the rescaled curves 
for a(t) according to Eq. (fl5)) with the extracted val- 
ues of the parameters are presented. As can be seen, 
the rescaled data generate a unified master curve, well 
described by the scaling form [Eq. (|15p]. So, crystalliza- 
tion kinetics of the system is defined by the time scales 
responsible for crystal nucleation and cluster growth, re- 
spectively. This has a similarity with results of Cavagna 
et al. for a lattice spin system [4(| [4l[ , where at low tem- 
peratures the fast nucleation of small and stable crystal 
droplets followed by slow activated crystal growth. 

Fig.[7]shows the 7-dependence of the nucleation rate i s 
at fixed temperatures. In this figure, the data obtained 
from the evolution of crystallization kinetics within 
Eq. (fT5|) are compared with the results of Yatsuoko- 
Matsumoto method at the known critical cluster size. 
Remarkably, both methods yield a very close behav- 
ior and reveal the same features for all the considered 
cases. As can be seen, i s (7) at the constant tempera- 
ture is nonmonotonic. Namely, the nucleation rate i s 
increases linearly at low strain rates. Then, at the val- 
ues 7 = 0.2r _1 -f- 0.8r _1 the nucleation rate levels off 
and reaches a maximum. On further increase of 7 nu- 
cleation rate starts to decrease. Interestingly, a similar 
7-dependence was earlier found by us for the phase trans- 
formation rate of the system [l2j . Results similar to those 
shown in Fig. [7] have also been observed for steady-state 
nucleation rate measurements against shear rate in very 
different systems, including an industrial polydisperse 
isotactic pol ypr opylene melt [20l | and a two-dimensional 
Ising model |16| . 

To understand the observed nonmonotonic behavior of 
is (7) we also consider the values of the product c g p s Ql 
and of the lag-time t c presented in Table Q] [the inverse 
lag-time vs. the strain rate is also shown in the inset of 
Fig. ©]. As can be seen, both quantities correlate di- 
rectly with nucleation rate, with a very similar nonmono- 
tonic variation. The rise of nucleation rate is therefore 
accompanied by the increase of the cluster growth and 
by the reduction of the time scale for transient regime 
and, vice versa, the decrease of nucleation rate occurs 
at retarding the cluster growth and increasing the lag- 
time scale. Such a correlation indicates directly that the 
changes in the nucleation rate are essentially a kinetic, 



rather than thermodynamic effect. 

Recently, the similar nonmonotonic behavior of nucle- 
ation rate vs. strain rate for a driven two-dimensional 
Ising model was revealed in Ref. [l6[. The authors had 
related the observed behavior with an interplay between 
shear-enhanced cluster growth, cluster coalescence and 
cluster breakup. It was additionally found in Ref. [l6| 
that shear-enhanced cluster coalescence and monomer at- 
tachment (single spin flip growth in Ising model) give the 
similar impact in the total ordering. Results of our study 
indicate rather that cluster growth and nucleation pro- 
cesses are merely correlated and their features are defined 
by kinetics of the transition. At the same time, a weak 
7-dependence observed for the critical cluster size [see 
Table |T] reflects the influence of shear drive on the par- 
ticle cohesion in a crystalline nuclei due to mechanical 
stresses. 



V. CONCLUSION 

In summary, we have performed nonequilibrium molec- 
ular dynamics simulations to study crystal nucleation 
and growth processes induced by shear drive in a metal- 
lic glass for a range of temperatures and strain rates. 
By applying a mean first passage time analysis, we de- 
fine the size of a critical cluster and the time scale of 
its appearance. We find that the nucleation-growth pro- 
cess has a non-stationary character, and the crystalline 
cluster grows with a time- and size-dependent rate. 

To describe the crystallization kinetics under shear ob- 
served in our simulations, the extension of the KJMA 
theory is suggested and compared with simulation data. 
As a result, an excellent agreement is obtained for all 
the considered cases. Further, we find that data for time 
evolution of the crystalline fraction at particular values of 
strain rates can be rescaled within this theoretical model 
to give an unified master curve. 

The observed nonmonotonic behavior of the nucleation 
rate J on the strain rate 7 at the fixed temperatures is 
very similar with that was recently reported for the case 
of a two-dimensional Ising model under shear [l6j . This 
behavior of J indicates directly that shear drive can speed 
up as well as suppress nucleation in a glass. Finally, we 
find that nucleation rate, lag-time and cluster growth are 
affected by the shear in a very similar and correlated way. 
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